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Abstract 

We summarize existing constraints on the epoch of reionization and discuss the observational probes 
that are sensitive to the process. We focus on the role large scale polarization can play. Polarization 
probes the integrated optical depth across the entire epoch of reionization. Future missions such 
as Planck and CMBPol will greatly enhance our knowledge of the reionization history, allowing us 
to measure the time evolution of the ionization fraction. As large scale polarization probes high 
redshift activity, it can best constrain models where the Universe was fully or partially ionized at 
early times. In fact, large scale polarization could be our only probe of the highest redshifts. 
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1 The Epoch of Reionization: How the Universe emerges 
from its dark ages 

1.1 Introduction 

Perhaps the most important goal of observational cosmology is to understand how structure arose 
in our Universe, how the stars, galaxies, and quasars around us today and at high redshifts formed. 
We currently believe that small density perturbations were imprinted onto the Universe during the 
inflationary era which then grew through gravitational instability, producing halos as well as the 
cosmic web of sheets and filaments. Measurements of the cosmic microwave background (CMB) 
anisotropics have been instrumental in establishing this picture (e.g. [if). 

To understand the CMB measurements or the large scale clustering properties of galaxies one 
only needs to understand the structure formation process in the linear regime. However, in order to 
understand how baryons collapsed into bound objects such as galaxies and galaxy clusters, and to 
determine how these objects affect their surroundings we need to take structure formation beyond 
the well-understood linear regime. At low or moderate redshifts {z < 6), galaxies and quasars can 
be observed directly with existing technology. However, the first generations of protogalaxies are 
not yet similarly accessible and the full population of galaxies and quasars at intermediate redshifts 
{z ~ 3) remains unknown. 

One approach to investigating the properties of galaxies and quasars is through their impact on 
the intergalactic medium (IGM). In particular, it is believed that radiation from these objects reion- 
ized the IGM and that this process depends on the distribution and nature of the ionizing sources. 
Thus, an understanding of reionization coupled with detailed measurements has the potential to 
significantly constrain models for the formation and evolution of galaxies and quasars. 

Observational evidence indicates that the IGM was reionized in two distinct phases. The first 
phase, completed by z ~ 6, led to the ionization of hydrogen as well as singly ionizing helium. 
However, the early sources of ionizing radiation would not have been energetic enough to doubly 
ionize helium (if they were stellar). The second phase of the reionization of the IGM, corresponding 
to the conversion of Hell into Helll, thus likely occurred later, when a sufficient population of hard 
sources (i.e. quasars) developed, and observations indicate that this happened at z « 3. 

Reionization encodes a wealth of information about the ionizing sources. For example, in the case 
of hydrogen reionization, the first sources of radiation likely formed in overdense regions, ionizing 
pockets of nearby gas. As more sources formed these HII regions grew and eventually overlapped. 
Because the structure of these HII regions and their growth depend on the nature, clustering, and 
abundance of the sources, the timing, morphology, and duration of reionization are highly sensitive 
to the pro perties of the first luminous objects and their relationship to the IGM (e.g.. Baas a 

saaH)- 



Empirical studies of the reionization of hydrogen (sometimes referred to as the Epoch of Reion- 
ization [EoR]) have placed constraints on this event but have not yet mapped it in detail. One 
constraint comes from the effects of the ionized gas on the CMB. The free electrons Thomson scat- 
ter CMB photons, washing out the intrinsic anisotropics but generating a polarization signal on large 
angular scales. The total scattering optical depth Tes is proportional to the column density of ionized 
hydrogen, so it provides an integral constraint on the reionization history. Recently, the Wilkinson 
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Microwave Anisotropy Probqj (WMAP) collaboration used five years worth of data to improve their 
determination of t^s • The combined analysis of WMAP and large scale structure data such as that 
coming from the Sloan Digital Sky Surve}o (SDSS) indicates that reionization probably happened 
in the redshift range Zr = 11 HA. More detailed information on the reionization history could be 
obtained by precisely measuring the large angular scale polarization [ll|, [l^, 13 1 . It is the goal of 
this white paper to discuss this possibility further. 

Along with the large scale CMB polarization measurements, we anticipate a wealth of other 
upcoming observations that will increase our knowledge of the EoR. The EoR and the preceding 'dark 
ages' - i.e., the period after cosmological recombination but before the birth of the first luminous 
sources - are the only remaining phases of structure formation as yet to be directly observed. As 
such, the study of the EoR is a research frontier for both observational and theoretical cosmology, 
and will grow as an increasingly active and exciting research field in upcoming years. We describe 
current constraints and future prospects from other observations in ^ and consider the way in 
which large scale CMB polarization measurements may complement these observations throughout 
this white paper. 



1.2 Motivating questions and how to answer them 

We believe that the study of reionization will shed light on the early process of structure formation. 
In practice we want to answer several specific questions about this early time in the history of 
our Universe. We want to determine which sources were responsible for reionization. Were they a 
population of early galaxies similar to present day ones? Were Population III stars important for this 
process? Did quasars or mini-quasars play a role? Did the sources of ionizing radiation sit in high or 
low mass galaxies? Did the ionizing radiation produced by the first sources have a strong feedback 
effect on the galaxy formation process, drastically changing the way galaxy formation proceeded 
thereafter? 

In some sense the study of reionization is akin to using the IGM as a laboratory for learning about 
the first sources and the earliest stages of the structure formation process. In order to understand 
how the study of the IGM can help us answer these questions, it is useful to think of two separate 
characterizations of the process. First we can think of the evolution of the average ionization fraction 
of hydrogen with time and we can subsequently consider the spatial fluctuations around that mean 
evolution. The latter is the study of the sizes of the ionized bubbles and the morphology of the 
ionized regions. This split is useful because different observational probes have different sensitivities 
to these different aspects. Furthermore the study of the mean and the spatial fluctuations are 
sensitive to different properties of the first sources. 

To learn about what type of halos the sources lived in one can study how rapid the evolution 
of the mean ionization fraction is. Sources that live in rarer more massive halos usually lead to a 
faster evolution of the mean. The time and duration of reionization also depends on the dumpiness 
of the IGM, as well as on radiative, chemical, and mechanical feedback from the first sources, which 
impact the formation of sources forming later on. An example of the type of information contained 
in the redshift evolution of the mean ionization fraction is shown in Figure 1, which demonstrates 

^See http://map.gsfc.nasa.gov/. 
^See http:/ /www. sdss.org/. 
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Figure 1: Ionization history of the IGM from p^]. The curves show example theoretical models 
describing the volume-filling fraction of ionized g function of redshift. Obser- 

vational probes of reionization aim to map-out, or constrain, the detailed dependence of 
ionization fraction with redshift. When does essentially the entire volume of the IGM 
become filled with ionized gas (xj = 1)? How extended is the reionization process? 
Answering these questions will provide significant constraints on the first generation of 
ionizing sources, and on the surrounding IGM. The example curves illustrate the impact 
of one relevant effect, radiative, thermal feedback, which acts to suppress the forma- 
tion of galaxies within ionized regions. The red solid curve ignores feedback, the black 
dashed line incorporates a 'maximal' level of thermal feedback, while the other curves 
include intermediate levels of feedback. (See [l4|] for details.) The timing and duration 
of reionization depend strongly on feedback processes such as this. 



the impact of radiative feedback. 

Furthermore any evidence of ionized hydrogen at very early times could signal the existence of 
an early and very efficient mode of star formation or perhaps an entirely new source of ionizing 
photons such as those that could come from annihilating or decaying dark matter. The study of 
large scale CMB polarization is probably the best and only probe of these very early stages of the 
ionization history and thus it is probably the only way to constrain some of these scenarios. 

The study of the spatial variations of the ionization fraction during reionization, that is deter- 
mining the size distribution of the ionized regions, also encodes a wealth of information about the 
sources responsible for reionization. A simulated model, illustrating spatial variations in the ion- 
ization fraction, is shown in Figure [21 The size distribution is mainly dependent on the clustering 
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Figure 2: Ionized bubbles during reionization from (15l] . Each panel is a 0.25 (co-moving) Mpc//i 
deep slice from a simulation box that is 65.6 Mpc//i on a side. The black regions show 
ionized bubbles and the color scheme is proportional to neutral hydrogen density. The 
panels show different stages of reionization, with Xi = 13%, 35% and 55% of the volume of 
the IGM ionized from top to bottom. Reionization studies aim to map-out, or constrain, 
the size distribution of these ionized bubbles at different stages of the reionization process. 
The left-most panels show results from a radiative transfer simulation of reionization, 
while the center and right-most panels show results from semi-analytic calculations. The 
blue points in the center and left panels show ionizing sources from the simulation. 

properties of the sources which in turn can help determine what mass halos those sources typically 
inhabit. Furthermore, and especially during the later stages of reionization, the presence of high 
density regions that act as photon sinks can set a maximum size for the ionized bubbles. Thus there 
is the potential to learn about the small scale density structure of the IGM. 
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1.3 Outline of this White Paper 

In section [2] we will describe different probes that can be used to constrain the ionization history of 
the universe. In section [3] we focus on the information that can be extracted from the large scale 
polarization of the CMB. In section U] we comment on the impact uncertainties in the ionization 
history have on the determination of other cosmological parameters. Finally for completeness in 
section [5] we discuss the measurement of secondary anisotropies. 

2 Survey of observational probes 

As mentioned in the introduction, a confluence of data sets from a variety of different observational 
probes should soon provide significant advances in our understanding of the EoR. In order to best 
understand how CMBPol may help with the study of reionization, we briefly describe several fu- 
ture probes, examining in particular how large scale polarization measurements may complement 
constraints from these other exciting studies. As we will see, upcoming observations should provide 
detailed constraints on reionization activity at redshifts close to and a bit higher than z ~ 6 while 
large scale polarization, and CMBPol in particular, may be our best way to study reionization at 
higher redshifts, z > 10 — 12, in the near future. 

We describe here four current and upcoming probes of reionization: high redshift quasar spectra, 
surveys for Ly-a emitting galaxies, optical afterglow spectra of gamma-ray bursts, and 21 cm ob- 
servations of the high redshift IGM. In ^we additionally discuss constraints from small scale CMB 
measurements, which will be made by several upcoming experiments, and potentially by CMBPol, 
provided its angular resolution is high enough. We also discuss the related topic of metal absorption. 

2.1 High Redshift Quasar Spectra 

Presently, our most detailed probe of the high redshift (z > 6) IGM comes from Ly-a forest ab- 
sorption spectra towards high redshift quasars. Regions with relatively large HI densities should 
appear as absorption troughs in quasar spectra, which presumably deepen and come to dominate as 
we approach the reionization epoch. Indeed, spectra of z ~ 6 quasars selected from the SDSS show 
some extended regions of zero transmission [16j. Unfortunately, the Ly-a cross section is very large 
and a relatively low average neutral fraction at the level of (-'^^hi) ~ 10~^ — 10~^ is sufficient to give 
complete absorption in the z ~ 6 Ly-a forest (l7| . This makes it difficult to distinguish a largely 
neutral IGM from a significantly ionized one with this probe. 

Nevertheless, several detailed studies have examined whether these spectra may indeed probe the 
IGM before reionization completes. These analyses have considered a range of different statistics: 
the redshift evolution of the average absorption in the Ly-a and Ly-/3 forest, the amount of scatter 
in the absorption from sightline to sightline, the abundance and extent of dark gaps in the high 
redshift Ly-a forest, the properties of the proximity zone regions in high redshift quasar spectra, the 
redshift evolution of metal abundances, and the temperature of the IGM from the z ~ 2 — 4 Ly-a 
forest. 

A plot of the redshift evolution of the Ly-a opacity is shown in Figure [3l More precisely, we 
show the redshift evolution of the effective Ly-a opacity, defined by r^g = —ln{F), where (F) is 
the average flux transmitted through the Ly-a forest. One can see that the opacity at z ~ 6 is 
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Figure 3: Redshift evolution of tlie Ly-a forest opacity from [l9|. The black points are measure- 
ments of the z < 5 Ly-a opacity from 18|]. The solid color ed p oints are measurements 
from individual Az = 0.15 stretches of Ly-a forest spectra [17|, with arrows indicating 
lower limits on stretches consistent with complete absorption. The black dashed and 
solid blue curves indicate different extrapolations of the z < 5.5 opacity to z >6. 



typically very large, t^s > 5, with some sightlines having only lower limits on the opacity - i.e., they 
are consistent with complete absorption in the 2; ~ 6 Ly-a forest. As mentioned before, complete 
absorption in the Ly-a forest at 2: ~ 6 only implies a weak constraint on the average neutral fraction. 
A slightly stronger constraint comes from spectra with complete absorption in the Ly-/? region of the 
forest. Since the Ly-/3 cross section is a factor of ~ 6 times smaller than that for Ly-a absorption, 
complete Ly-/3 absorption implies a more stringent lower limit on the neutral fraction than complete 
absorption in Ly-a, but the resulting constraints are still consistent with a mostly ionized IGM [17| . 

A second feature of this figure, however, is that the opacity appears to evolve rather sharply with 
redshift around 2 ~ 6. Perhaps the sharp evolution in the opacity is signaling correspondingly rapid 
changes in the ionization state of the IGM? This is currently a matter of debate. In 171] the authors 
fit a power-law to z < 5.5 data (shown as the black dashed line in Figure [3]), and find that the 
z > 5.7 data show a rapid upward departure from this fit fl^. Using a simulated model for the IGM 
density distribution [2d|], they argue that this upward departure requires very rapid evolution in the 
ionization state of the IGM. On the other hand, the authors of |19( | use an empirically motivated fit 
to the opacity distribution and claim that rapid changes in the ionization state are not required 19|] . 
The blue solid line in Figure [3] shows an extrapolation of their model from z < 5.4 data to higher 
redshift, assuming a slowly varying ionization state. This model also fits the z > 6 Ly-a opacity 
data. The main issue here is that, even in a highly ionized IGM, the Ly-a absorption cross section is 
so large that only quite underdense regions allow transmission through the 2; ~ 6 Ly-a forest. This 
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means that the opacity evolution is sensitive to the precise abundance of underdense regions, as are 
the constraints on the ionization state ( 2l| , 2^ ) . Further theoretical modeling may help resolve 
this debate. 

In addition to the large average absorption, the scatter in the absorption from sightline to 
sightline is quite large at z ~ 6 Indeed, the colored points in Figure E] show order unity variations 
in the opacity, after averaging over A2; = 0.15 or ~ 50 co-moving Mpc/h. These variations might 
indicate incomplete and patchy reionization with some sightlines piercing through mostly ionized 
gas, and others passing through long neutral stretches 23|]. However, variations in the line of sight 
density field give rise to surprisingly large opacity fluctuations at 2: ~ 6, and so it is difficult to 
discern, on the basis of this scatter, whether these spectra really probe the IGM before reionization 
completes 22]. Current measurements appear at least broadly consistent with density fluctuations 
alone 0,[2|]. 



The tightest constraints claimed on the ionization state of the 2: ~ 6 IGM come from measure- 
ments of the proximity regions around 2; ~ 6 quasars. Several authors, starting with 25|], have 
argued that these re gion s are small, indicating that quasar ionization fronts are expanding into a 
largely neutral IGM 25|, |26|, |27|. However, subsequent more detailed studies have shown that it is 
hard to distinguish the case of a quasar ionization front expanding into a partly neutral IGM, and 
the case of a mere proximity zone in a mostly ionized IGM (e.g., [28l . |29||). Detailed modeling of 
inhomogeneous reionization, incorporating the tendency for quasars to live in overdense regions that 
are ionized before typical regions, is required to best interpret these observations 291 ]. 

One can also measure the redshift evolution of metal line abundances in high redshift quasar 
spectra, and use this to constrain the reionization history. The mass density in Carbon IV appears 
to be a surprisingly flat function of redshift from 2 ~ 2 — 6 30, l3ll] , requiring significant levels of star 
formation and metal enrichment prior to 2; = 6. Another metal line tracer is the 01 line which has an 
ionization potential similar to that of hydrogen, lies redward of Ly-a, and has a significantly lower 
optical depth than Ly-a, making it a potential probe of the IGM before reionization completes 32l |. 
In fact, high-resolution Keck spectra of the 2 ~ 6 SDSS quasars do reveal some 01 lines, with 4 out 
of 6 detected systems lying towards the highest redshift quasar presently known 33|]. Interestingly, 
some of the 01 systems are nearby regions that show transmission in the Ly-a and Ly-/3 forests of 
this quasar [33]. The interpretation of these observations is unclear: the 01 systems might reflect 
dense clumps of neutral gas in a highly ionized IGM, or instead could indicate inhomogeneous metal 
pollution in a more neutral IGM. 

A final constraint comes from measurements of the temperature of the Lyman-a forest at 2 ~ 2- 
4, which suggest an order unity change in the ionized fraction at Zr < 10 
argument depends on the timing and history of Hell reionization (e.g., 36l. 137 

Some progress will be made in the near future by detecting more quasars, some at higher redshifts 
(2 > 6.5), using widefield, deep near-IR surveys, although the rarity of bright background quasars 
at these redshifts makes this challenging. Improvements in modeling inhomogeneous reionization 




35(, although this 



should also help our understanding of 2 > 6 absorption spectra [29|, |38|, |39[ . However, the large Ly-a 
opacity expected in a mostly ionized medium is clearly a fundamental limitation, and it is likely 
impossible to study all but the final stages of reionization with this approach. 
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2.2 Lyman-alpha Emitters 



An effective way to find high redshift galaxies is to search for their Ly-a emission hnes, since young 
galaxies frequently have strong Ly-a emission 40|]. Ly-a emission searches have an advantage over 
high redshift Lyman break surveys in that they target narrow wavelength intervals, in between strong 
night sky background lines, in search of prominent emission lines. This allows one to detect galaxies 
that are unobservable by Lyman break selection owing to the strong night sky background at relevant 
wavelengths. Existin g su rveys have discovered several hundred galaxies with this techni que (e.^ 



41 



42 



43 



44, 45, 46, 47, 48, 49]), and other programs will begin taking data soon (e.g., [sol-lsil. Is^]). 
The largest existing hi gh r edshift sample from the Subaru Deep Field consists of approximately 50 
galaxies above z 



6.5 



X:= .3 



X; = .5 



x-= .7 




Figure 4: Ly-a emitting galaxies during the EoR. The top panel shows maps of the ionization field 
from a simulation at three different (volume-averaged) ionized fractions: (xj) = 0.3, 0.5, 
and 0.7 from left to right respectively. The middle panel shows the true distribution 
of Ly-a emitting galaxies from the corresponding simulation slices. The bottom panel 
shows the apparent distribution of the same galaxies, based on observing the sources' 
Ly-a emission lines. From 53|. 
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During the EoR, neutral hydrogen gas in the IGM will extinguish galactic Ly-a emission lines. 
Specifically, the optical depth to Ly-a absorption in a significantly neutral IGM is so large that even 
the red side of the Ly-a line is attenuated by damping wing absorption [5j]. As a result, galaxies that 
lie sufficiently close to the edge of an HII region, next to a significant column of neutral hydrogen, 
will be unobservable in Ly-a. On the other hand, galaxies which reside towards the center of a 
sufficiently large HII region avoid complete attenuation, and will still be visible in Ly-a emission. 
Since large HII regions form around high density peaks, one expects the clustering of Ly-a emitters 
to increase dramatically towards the start of reionization 53l. I55l.l56|. In conjunction, the abundance 
of Ly-a emitters should drop off rapidly as one looks back deeper into the EoR. 

These trends can be seen in the simulation slices of Figured! Comparing the top, middle, and 
bottom panels, one can see how the presence of ionized bubbles modulates the abundance of emitters 
observable in Ly-a: sources close to the center of large ionized regions are observable, while other 
sources are attenuated out of the Ly-a sample. Comparing the left-most and right-most panel, it 
is clear that the suppression is much more significant in the earliest stages (here when the volume 
ionized fraction is Xi = 0.3), compared to later stages (here when Xi = 0.7). It is also clear that 
- since the ionized regions quickly grow rather large - one needs to observe a large field of view 
to reliably measure the abundance and clustering of Ly-a emitters during reionization. To best 
constrain reionization with this technique, one hence wants a widefield survey that is deep enough 
to find large samples of emitters in the early stages of reionization. 

Several authors have, in fact, used existing Ly-a emitter samples to place constraints on the 
ionization state of the IGM. Malhotra & Rhoads 57| constrained the ionized volume fraction in 



the IGM by counting the abundance of Ly-a emitters at z = 6.6, and requiring a minimum ionized 
volume around the sources to avoid attenuating their Ly-a photons. From this argument, they 



find that more than 50% of the volume of the IGM is ionized at z = 6.6. Kashikawa et al. [46 1 
observe a ~ 2 — cr decline in the bright end of the luminosity function between z = 5.7 and z = 6.5, 
which they attribute to incomplete reionization at z > 6. Several other authors argue, however, 
that the observed evolution is consistent with evolution in a post-reionization IGM [ssl . 59]. 



Finally, McQuinn et al. [53t | use the clustering of the z = 6.6 Subaru Deep Field sample to exclude 
a significantly neutral Universe with Xi < 0.5 at 2 — cr confidence. 

Further progress can be expected, with numerous planned surveys attempting to find Ly-a 
emitters at high redshift. In particular, an extension to the existing Subaru Deep Field survey 
should extend measurements to z = 7.3 over a ~ 2 deg^ field of view. This amounts to a factor of 
~ 10 boost in field of view compared to the existing Subaru Deep Field. It will likely be difficult 
to push this method to significantly higher redshifts, since strong night sky emission lines force one 
to make near IR measurements from space, and since we require observations over a wide field of 
view. We expect this method to be most interesting in the near future if a significant fraction of the 
volume of the IGM is still neutral near z = 7. 



2.3 GRB Optical Afterglows 



Another possibility is to search for damping wing absorption in high redshift GRB optical afterglow 



spectra [54, ISO, l6l|) These sources are extremely luminous, can potentially be detected at 



very high redshift, and have the important advantage of a simple power-law intrinsic spectrum. 
Additionally, GRBs are more likely to exist in typical HII regions than quasars or the massive. 
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highly luminous galaxies that are easiest to detect at high redshift. Thus it may be possible to 
identify damping wing absorption in individual GRB afterglow spectra. 
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Figure 5: GRB optical afterglow spectra and DLAs. The blue line shows a portion of the observed 
spectrum of the GRB optical afterglow spectrum of GRB050904, at z = 6.3 (from 611]). 
The black solid line is a model spectrum, consisting entirely of absorption from a DLA 



associated with the GRB host galaxy, of column density A'^hi = 10 



21.6 



cm 



The red 



dashed line is a model that includes both DLA absorption and damping wing absorption 
from a completely neutral IGM. The presence of a large column density local absorber 
makes it extremely difficult to tell whether there is absorption from neutral hydrogen in 
the surrounding IGM. From 62] • 



A difficulty with this probe of reionization, however, is that many optical afterglows show damped 
Ly-a absorption from neutral hydrogen in the GRB host galaxy, at a level that is large enough to 
overwhelm the damping wing opacity from neutral gas in the surrounding IGM. For example, the 
highest redshift GRB afterglow detected thus far, at z 



6.3 



m 



shows evidence for a strong 
Damped Ly-a Absorber (DLA) in its host, with a column density of Ahi = 10^^'^cm^^. This 
absorption (shown in Figure [5|) is so strong that it is difficult to tell whether or not there is damping 
wing absorption from the surrounding IGM. The damping wing profile from extended neutral gas in 
the IGM has a different wavelength dependence than that from a compact, local absorber, and we 
can therefore still hope to unambiguously detect neutral gas in the IGM using GRB afterglow spectra 



541 ]. To detect the IGM damping wing, a GRB with a weaker local absorber, smaller than roughly 



A^Hi ^ 10 cm~ , is needed. In addition, a high signal-to-noise spectrum is desirable, demanding 
rapid near IR follow-up of probable high redshift bursts. Some afterglow spectra at lower redshift 
show local absorbers with significantly smaller column density than the z = 6.3 burst, and occasional 
spectra show no DLAs whatsoever [63(. We may therefore hope to find a clean line of sight towards 
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a GRB at higher redshift and apply this test. 

While GRB optical afterglow spectra may provide a good way of detecting the presence of 
neutral hydrogen in the IGM, perhaps even at very high redshift, it will be difficult to extract the 
ionized fraction from such measurements. The challenge here is that reionization is expected to 
be extremely inhomogeneous and the level of damping wing absorption should vary considerably 
towards one GRB host or another. Some GRBs will reside in large HII regions even rather early 
in the reionization process, since large ionized regions can form under the collective influence of 
highly clustered neighboring galaxies. The GRBs in large HII regions will suffer little damping 
wing absorption, while GRBs in smaller ionized bubbles or close to the edge of a large bubble will 
suffer more damping wing absorption jo^] . This is illustrated in Figure [U which shows that even 
when the IGM is as much as half neutral, there will be significant scatter in the level of damping 
wing absorption. This is similar to the effect depicted in Figure [U but here we are attempting to 
extract information based on the spectrum of an individual object, rather than statistical trends 
based on the abundance and clustering of an entire population of galaxies. Extracting the ionized 
fraction from this method may require a prohibitively large sample of high redshift afterglows. We 
are hopeful that this method will definitively detect neutral gas at high redshift, indicating that 
reionization is not yet complete at a given redshift, but expect other probes to best map out the 
redshift evolution of the ionized fraction. 



2.4 The 21 cm Line 



One of the best ways to study reionization may be to detect redshifted 21 cm emission from neutral 
hydrogen gas in the high redshift IGM. Whenever the spin temperature of neutral hydrogen is 
different from the CMB temperature, neutral hydrogen atoms can be seen either in emission or 
absorption against the CMB at the redshifted wavelength of their 21 cm transition. After the first 
sources appear, the gas is expected to have been heated and the kinetic temperature of the gas 



should exceed that of the CMB 6J]. Furthermore, Ly-a photons from the first sources should 



couple the spin temperature of the 21 cm transition to the gas temperature long before the IGM 
is significantly ionized [65[. Thus, hydrogen should appear in emission during the EoR and high 
resolution observations of the 21 cm transition as a function of both frequency and angle can provide 
a three-dimensional map of reionization (e.g., 0, 67]). Several experiments are now underway or 



being planned to measure this signal, including the 21 cm ARRAY to be operated in China (formerly 
known as PAStH), the Low Frequency Array (lofaeS), the Murchison Widefield Array (MWA@), 
the PAPERr experiment, an effort at the GMRlQ telescope, and ultimately the Square Kilometer 
Array (SKAo). This method has several key advantages: the optical depth is low and so one does not 
suffer from the saturation problems that plague Ly-a forest studies, it does not require a high redshift 
background source, and the probe is a spectral line, allowing one to extract full three-dimensional 
information about reionization. For a review, see 68 1. 



•^http:/ /web.phys.cmu.cdu/^past/ 
^http:/ /www. lofar.org 

^ http://web.haystack.mit.edu/MWA/MWA.html 

^http:/ /astro. berkeley.edu/ dbacker/EoR/ 

^Ue-Li Pen private communication. 

® http://www.skatelescope.org for details on the SKA 
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Figure 6: Sample variance and GRB afterglow spectra from [64]. Top panel: The lines show the 
neutral fraction along several lines of sight towards a GRB host, when the IGM is ~ 50% 
ionized by volume. Bottom panel: The damping wing absorption vs. wavelength for 
each of the nine sightlines shown in the top panel. Some lines of sight lead to significant 
levels of damping wing absorption, while other lines of sight result in very little damping 
wing absorption. It is hence difficult to discern the ionized fraction from a few optical 
afterglow spectra, even if they unambiguously show damping wing absorption from the 
IGM 
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Futuristic experiments like the SKA should have the sensitivity to produce relatively detailed 
maps of the reionization process. These maps, taken at several different frequencies, will amount 
to a reionization 'movie': they will depict the growth of HII regions around individual sources, 
their subsequent mergers with neighboring HII regions, and detail the completion of the reionization 
process, whereby the entire volume of the Universe becomes filled with ionized gas. Reionization 
maps from the SKA should rather directly reveal the size distribution of ionized regions during 
reionization, and the redshift evolution of the ionized fraction. First generation 21 cm surveys, 
such as the MWA, LOFAR, GMRT, and the 21 cm ARRAY, will not have the sensitivity to make 
detailed maps of the reionization process, but will allow for a statistical detection j69f|. For example, 
the MWA should measure the power spectrum of 21 cm fluctuations over roughly a decade in scale, 
in each of several independent redshift bins 70(]. These measurements will already be quite valuable, 
but may not be entirely straightforward to interpret. 




Figure 7: Sensitivity of the MWA for measuring the 21 cm power spectrum. The red curves show 
dimensionless 21 cm power spectra from a reionization model at different redshifts and 
ionization fractions, as labeled. The red error bars indicate statistical error estimates for 
1, 000 hours of observations. The blue error bars indicate error estimates for a different 
conflguration of the MWA antennas. Foreground cleaning will likely prohibit power 
spectrum measurements at wavenumbers smaller than indicated by the black dashed 
lines. From [to! ]. 
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Here we briefly discuss what we might hope to learn from these first generation experiments, 
which will come online before CMBPol. A model for the 21 cm power spectrum during reionization 
is shown in Figure [71 along with approximate error bars for one year of MWA observations 70(] . 
At early times - as seen in the top left hand panel of the figure - before significant ionized regions 
grow in the model, the 21 cm power spectrum traces the density power spectrumjf] As large ionized 
bubbles grow around early generations of galaxies, the amplitude of the 21 cm power spectrum 
increases on large scales, and the slope of the 21 cm power spectrum on MWA scales flattens. 
This behavior is seen in the bottom left hand panel of Figure [71 Finally, the amplitude of 21 cm 
fluctuations drops as the Universe becomes still more ionized and neutral hydrogen becomes scarce 
(bottom right hand panel). The amplitude on MWA scales is typically maximal around the time 
when the Universe is half ionized. To recap, one expects the amplitude of the power spectrum to rise 
and fall as the Universe becomes ionized, and the slope of the power spectrum to flatten as ionized 
regions become large and impact the scales probed by the MWA. These trends should be robust, 
but detailed modeling is clearly required to best extract information about the redshift evolution of 
the ionized fraction and the size distribution of ionized regions during reionization. Note that the 21 
cm signal should be highly non-Gaussian during the EoR: there is more information than contained 
in the power spectrum alone and there may be better ways to extract constraints on the ionized 
fraction and bubble sizes from the observations. 

The flgure also shows error bar estimates for about one year of observations with the MWA at 
different redshifts. On small scales, the detector noise becomes large because the instrument has few 
antennas at long baselines and it's high k sampling is thus mainly limited to modes in the frequency 
direction. Large scale modes, on the other hand, will be lost to foreground cleaning. Astrophysical 
foregrounds are a considerable challenge for these observations, as the foregrounds are expected to 
be four orders of magnitude larger than the signal itself. Fortunately, known astrophysical fore- 
grounds are spectrally smooth and hence distinguishable from the high redshift 21 cm signal itself 
(e.g., 72]). Nonetheless, one inevitably loses long wavelength modes along the line of sight in the 
foreground cleaning process. A rough estimate is that foregrounds preclude measuring 21 cm power 
at wavenumbers short of the dashed line shown in the figure. This implies that the MWA is sensitive 
to roughly a decade in scale, probing wavenumbers between A; ~ 0.1 — l/i Mpc~^. The red and blue 
points show error estimates for different configurations of the MWA's antennas (see [70| for details). 

Another important point for our present discussion is that the sensitivity is a very strong function 
of redshift, likely prohibiting upcoming measurements from probing very high redshifts. The reason 
for this is that the radio sky is much brighter at low frequency, with the sky brightness scaling as 
Tsky oc (1 + z)^'^, which results in a rapid increase in detector noise towards low frequencies, and in 
limited high redshift sensitivity. This trend is evident in Figure [71 Part of the reduced sensitivity 
towards high redshift results because reionization occurs relatively late in this model, and the signal 
is weak at high redshift here. The signal is typically maximal in the middle of reionization, and the 
high redshift detectability will be boosted in models where the middle of reionization occurs at higher 
redshifts than considered here. The increased noise at high redshift is unavoidable, however, and 
upcoming surveys will not be sensitive to very high redshift stages of reionization, losing sensitivity 
at z > 10 — 12. Moreover, because of the loss of sensitivity towards high redshift, first generation 



^Spin temperature fluctuations are neglected in this model. Including these will modify the predictions in 
the early stages of reionization (e.g., [7l|). 
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surveys will likely focus their efforts exclusively on relatively low redshifts. The enormous data rates 
involved with correlating large numbers of antenna tiles will limit the observing bandwidth of the 
first generation experiments. The MWA, for example, will process an instantaneous bandwidth of 32 
Mhz, corresponding to a redshift extent of Az = 1.8 near z ~ 8. The full spectral range covered by 
the instrument is much larger, 80 — 300 Mhz, but one does not observe the full spectral range for free. 
This fact, and the reduced sensitivity towards high redshift, mean that initial MWA observations 
will likely focus on contiguous 32 Mhz stretches near z > 6. 

To summarize, we expect exciting constraints on reionization from several independent probes in 
the next few years. These probes should help determine the redshift evolution of the ionized fraction, 
and the size distribution of ionized regions during reionization, which will in turn constrain models 
for the ionizing sources and early structure formation. A common feature of each of these probes 
is that constraining reionization activity at redshifts close to z ~ 6 is much easier then probing 
reionization at higher redshifts, z > W — 12. 



3 Large Scale CMB Polarization 

The latest WMAP measurements of the large scale CMB polarization constrain the redshift of 
reionization to be Zreion = H-O ± 1.4 {731], assuming instantaneous reionization. As large scale 
polarization measurements improve, this assumption may be relaxed and further details of the 
ionization history will be constrained by the data. We explore how this may come about in the 
current section. We start by discussing constraints on the total optical depth and progressively 
include more parameters to assess the amount of information that a future mission such as CMBPol 
could obtain. 



3.1 Future Constraints on the optical depth 

Given the improved sensitivity, Planck is expected to be cosmic variance limited in polarization out 
to I ~ 10, improving on the determination of r from WMAP5 by a factor 2.5: a-^ = 4.7 x 10~^ 



7J]. This improvement does not significantly depend on whether a general reionization history or 



mm- 



an instantaneous reionization scenario is assumed 

Planck results may be improved with the help of sub-orbital experiments observing a large 
fraction of the sky with low sensitivity. Indeed, combining Planck with a cosmic-variance limited 
experiment out to / = 20 — 30 — 50 and covering half of the sky would further reduce the error on r to 
CTt- = 4.2 — 4.0 — 4.0 X 10^'^. In this case, sky coverage rather than sensitivity is the main limitation. 
More realistic calculations for planned experiments suggest that these expectations are not affected 
by actual noise properties or scanning strategies [ttJ. 

Assuming an instantaneous reionization process a full-sky cosmic variance limited experiment 
to / = 2500 on would have 0"t- = 1.9 x 10^^: more than a factor two lower than what Planck will 
obtain. Non-negligible noise and beam effects can degrade this performance; results for proposed 
configurations ([TS!]) are presented in table [H 
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Figure 8: Marginalized distributions for r and relat ed p arameters, for a fiducial model with instan- 
taneous reionization and r = 0.087 from [78]. Curves show the expected performance of 
Planck (green/solid lines), CMBPol (the EPIC 2m configuration, red/dashed lines) and 
a full sky cosmic variance limited experiment (yellow/dotted lines). For both Planck and 
CMBPol we assumed a 80% sky coverage. Curves include the effect of marginalization 
over the minimal ACDM parameter set. 



Planck 


CMBPol 


Cosmic Variance 


T 4.7 X 10-3 


2.5 X 10-3 


1.9 X 10-3 


As 9.4 X 10-3 


5.1 X 10-3 


4.3 X 10-3 


as 6.7 X 10-3 


3.3 X 10-3 


2.6 X 10-3 



Table 1: MCMC error estimates for r and related parameters, assuming instantaneous reionization 
78| 



3.1.1 Constraints on the timing of reionization 

The tight error bar on r expected for CMBPol should translate into a strong constraint on when 
reionization happens, constraining models for the ionization history and early structure formation, 
such as those of Figure 1. In order to illustrate this, we consider a simple two parameter model for 
the ionization history: the redshift of reionization, z^, is the redshift at which the hydrogen is half 
neutral, and "the duration of reionization," A, which is the width of a tanh function that describes 
its evolution: 



_ ^Of) = / ] 1 + tanh 



(l + ^^)3/2_ (1+^)3/2 



A 



(1) 



where / = 1 + nn^/nH ~ 1.08. This is currently the default parametrization of CAMB code. 
Note that there is no real physical significance to this parametrization, it is merely mathematically 
convenient. See Eq.(95) on pp.14 of CAMB noteJ^ for more details. Note that CAMB adds electrons 
from ionization of Hell, which increases Xe{z) by about 10% at z < 3.5. This has a negligible 
contribution to the optical depth, however, since most of the contribution to Tg comes from higher 
redshifts. 

We then examine how tight a constraint on reionization redshift, Zr (defined here as the redshift 
at which the IGM is half neutral), CMBPol can place, marginalizing over the uncertain duration of 
reionization, parametrized here by A. We contrast this constraint with one from existing WMAP 
5 year data. The result of this calculation is shown in Figure [9l for two different choices of input. 



^http:/ /cosmologist.info/notes/CAMB.pdf 
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Figure 9: Improved constraints on reionization liistory with CMBPol. Top left: Tlie solid line 
shows constraints on the reionization redshift from WMAP data, after marginalizing over 
the duration of reionization. The dashed lines are estimates for the constraints obtainable 
with CMBPol in each of a high redshift and a low redshift reionization model. Top right: 
The ionization history in each fiducial model. Bottom left: The TE cross power spectrum 
in each fiducial model, along with CMBPol error estimates. Bottom Right: The same for 
the EE power spectrum. 
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'true' reionization history. The two models are a low reionization redshift (zr = 7, A = 0.5) model, 
and a high reionization redshift {zr = 13, A = 0.5) model. Both models are consistent with current 
WMAP data to within a 99% confidence level. In each underlying model, the figure illustrates that 
CMBPol should provide a tight constraint on the timing of reionization. The tight constraint on 
high Zr models may be very valuable, if the IGM is indeed reionized at high redshift. As emphasized 
in ^ other observational probes are unlikely to be sensitive to high redshifts, z > 10 — 12, soon. 
In this case, others probes may first deduce that there are significant quantities of ionized gas at 
Zr ^ 10 (by failing to detect neutral gas at lower redshifts), but will learn little else. The figure 
illustrates that CMBPol can obtain much more detailed information in these high Zr scenarios. If the 
IGM is reionized at lower redshift, other observational probes may already have definitively detected 
diffuse neutral gas in the IGM. In this case, CMBPol measurements should still tighten constraints 
and provide an independent check with a different set of systematics. 



3.1.2 Constraints on high and low redshift contributions to the optical depth 

The improved sensitivity and resolution of CMBPol will push E-mode power spectrum measurements 
into the i ~ 10—20 range, and allow one to constrain more than just the total optical depth. The next 
best measured quantity will be the difference between the high redshift and low redshift contributions 
to the total r [79| (e.g., the shape of the second principal component (PC) of Xe{z) in ^3.1.5p . Given 
that the current constraint on instantaneous reionization models places the redshift of reionization 
at z « 10 ( 

^reion — 11.0 it 1.4 from WMAP5 [73| ), it is interesting to ask how well polarization 
can limit the contribution to the total optical depth from z < 10 and from z > 10. Instantaneous 
models maximize the low-redshift contribution, so extended reionization and more exotic models 
with partial ionization at high redshift would show up as detections of nonzero t{z > 10) and a 
reduced value of t{z < 10) relative to the total optical depth. 

The principal component analysis of § 13.1.51 provides one way to compute general bounds on the 
optical depth from fixed wide bins in redshift. The optical depth due to ionization between redshifts 
zi and Z2 is 

t(zi <z<Z2) = 0.0691(1 - Yp)nbh / dz I. .,' Xe{z). (2) 

Jzi H[z)/Hq 

By computing Xe{z) from the PC amplitudes of individual MCMC samples [Eq. ([5|)], we can use 
eq. ([2]) to obtain constraints on the derived parameters t(6 < z < 10) and t(10 < z < z^^x). [Since 
we assume full ionization at z < 6 for all models, the constant contribution of t{z < 6) ~ 0.036 is 
subtracted from the low redshift optical depth parameter.] 

For a given choice of Zmax we expect parameter constraints to only be accurate as long as there 
is no significant ionization between Zmax and z ~ 1100. Nevertheless, estimates of t(10 < z < z^ax) 
can be sensitive to ionization events at z > Zmax due to the fact that high redshift ionization affects 
a wide range of multipoles that overlap the range affected by ionization at z < Zmax- Therefore, 
r(10 < z < Zmax) can effectively be thought of as t{z > 10) [79(|. If future data show signs of high 
redshift ionization, then it will be necessary to check the robustness of constraints by increasing the 
chosen value of Zmax to ensure accurate parameter estimates. 

The constraints on the optical depth coming from low redshift, r(6 < z < 10), and from high 
redshift, t{z > 10), are plotted in Fig. [10] for current data and for the expected sensitivity of Planck 
and cosmic variance limited CMBPol. The error forecasts assume the same fiducial ionization 
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history and cosmology as for Fig. [T5| which has roughly equal optical depths of ~ 0.02 coming from 
6 < z < 10 and z > 10. 

The contours in the left hand panel of Fig. [10] are elongated in the direction of constant r, 
showing that the total optical depth is much better constrained than its contributions from more 
restricted redshift ranges even in the cosmic variance limit. The contours are cut off at the upper 
edge of the plot in Fig. [TU] due to the upper limit t(6 < z < 10) < 0.039, which comes from requiring 
that the ionized fraction be no more than unity over this range of redshifts. 

While the maximum likelihood models from the Planck and CMBPol MCMC analyses coincide 
with the fiducial values plotted as the cross in Fig. [TOl the confidence regions are shifted slightly in 
the direction of higher optical depth coming from z < 10. This may be partly due to the fact that 
6 < z < 10 is a fairly narrow redshift range compared to the full extent of the principal components, 
which we have taken to be 6 < z < 30 here. Consequently, some of the high variance PCs that can 
safely be neglected when considering Cf^ or total r due to having many oscillations in redshift are 
much smoother over the range 6 < z < 10. By fixing these components to have zero amplitude, 
we are ignoring any potential perturbations they can induce in r(6 < z < 10). The maximum 
perturbation to t(6 < z < 10) from the next principal component (p = 6) is ~ 0.005, comparable 
to the 68% uncertainty in this parameter with cosmic variance limited data. However, reionization 
histories with such large values of the higher-variance PCs are extreme and possibly even unphysical, 
so in practice the effects of the components that we ignore on the constraints in Fig. [10] should be 
reasonably small. 

The key point illustrated in Figure [TOl is that only CMBPol can strongly reject a significant high 
redshift contribution to the total optical depth. A strong limit on the high redshift contribution to 
the optical depth would rule out the presence of an early and efficient mode of star formation, and 
constrain the possible contribution of ionizing photons produced by annihilating or decaying dark 
matter. Alternatively, detecting a significant high redshift contribution would point to an exotic 
source of ionizing photons. 

3.1.3 Constraints on the duration of reionization 

Next we consider whether large scale polarization measurements can constrain the overall duration 
of the EoR, in addition to the redshift at which the IGM becomes half-neutral. 

We use the simple two parameter model described in section 13.1.11 We explore two sets of 
parameters: 

(i) "Low r model," (z^,A)=(10.5, 2.0), whose optical depth, r = 0.088, is close to the WMAP 
best-fitting value, twmap = 0.087 ± 0.017 173], 

(ii) "High T model," {zr,A)={13.0, 3.5), whose optical depth, r = 0.119, is close to the upper 95% 
CL limit of the WMAP value. 

Figures [11] and [12] show the WMAP 5-year constraints (blue) and forecasts (yellow for Planck 
and orange for CMBPol) for the low and high optical depth scenarios. In the low r model although 
Planck and CMBPol tightly constrains Zr (as in ^3.1.ip . the duration of reionization. A, is essentially 
unconstrained. Even the unphysical case of an instantaneous reionization model lies within the 1 — a 
error contour. Models with large A are either within the 2 — a error contours, or are already excluded 
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Figure 10: Constraints on contributions to the reionization optical depth from 6 < 2; < 10 and 



z > 10 with temperature and polarization data from WMAP5 {blue) and forecasted for 
large scale polarization from Planck (gold) and CMBPol (red). Left: 68% and 95% CL 
2D contours computed using the same set of MCMC samples as in Fig. [T5l The values 
of t(6 < z < 10) and t(z > 10) for the fiducial ionization history are indicated by a 
cross. Right: Marginalized posterior probability {solid) and mean likelihood {dotted) 
of t{z > 10). The discrepancy between the posterior distribution and mean likelihood 
for WMAP and for Planck at low t{z > 10) is due to the priors set by the range of 
possible values of the optical depth from 6 < 2; < 10 and z > 10. Only for CMBPol is 
there both a small number of MCMC samples and a small likelihood at t{z > 10) = 0. 



by z ~ 6 quasar spectra (see §2.1[) . This is partly because the ionization history in this model is 
symmetric around Zr- For instance, our current parametrization does not allow models where the 
ionization history rises abruptly for z < Zr, but is more gradual above z > Zr. In this case, one can 
constrain the high redshift contribution to the opacity, as we demonstrated in ^3.1. 2[ In our high r 
input model, one expects a weak constraint on the duration parameter. A, from CMBPol. In this 
case, the instantaneous reionization and very rapid reionization models can be rejected at 2 — o" by 
CMBPol, although Planck can not distinguish such models. Unfortunately, it is hard to distinguish 
between more probable models with, e.g., A = 2 and A = 4. 

3.1.4 Constraints on a physical model 

The 'zr- — A' parametrization of the previous subsection is convenient, but somewhat arbitrary, and 
so we would like to explore an additional model to test whether this impacts our main conclusions. 
In this subsection we consider an alternate two parameter model for the reionization history. This 
model is also somewhat simplistic, but in combination with our previous model and the principle 
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Figure 11: Low r model, Zr = 10.5 and A = 2.0. Planck or CMBPol cannot distinguish between 
the extended and instantaneous reionization scenarios for this case. {Top Left) Pro- 
jected constraints on the epoch and duration of reionization, parametrized by the tanh 
model given in Eq. ([1]). The blue contours show the existing WMAP 5-year constraints 
at the 68% and 95% CL, while the yellow and orange contours show the projected 
constraints from Planck and CMBPol, respectively. The diamond and star show the 
input model and an instantaneous reionization model with A = 0.5, respectively. {Top 
Right) Reionization history. The solid line shows the input model (diamond in the top- 
left panel), while the dashed line shows an instantaneous reionization model (star in 
the top-left panel). {Bottom Left) Temperature- E' polarization cross power spectrum. 
The boxes show the expected 68% CL bounds from CMBPol. The solid and dashed 
lines correspond to those in the top-right panel. {Bottom Right) E polarization power 
spectrum. The boxes show the expected 68% CL bounds from CMBPol. The solid and 
dashed lines correspond to those in the top-right panel. 
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Figure 12: Same as Fig. llll but for the high r model, z^. = 13.0 and A = 3.5. Although Planck can- 
not distinguish between the extended and instantaneous reionization scenarios, CMBPol 
can reject the instantaneous reionization at more than 95% CL. These models are dis- 
criminated best at Z ~ 20. 



component analysis of §3.1.51 it should provide a good sense for CMBPol's capabilities. 

In the model of this section, the ionization fraction is directly related to the fraction of mass 
that is in dark matter halos above some minimum mass. This form is motivated by assuming that 
the star formation rate is proportional to the rate at which matter collapses into galactic host halos. 
See e.g. 



68 



80| for a discussion of this model. More specifically, we assume that 



Xe{z) = 1 — exp — C erfc 



(3) 



Here the error function term, is the collapse fraction for halos with mass above Mmin according 
to Press-Schechter theory, 5c{z) = 1.686/Di(z) is the linear theory overdensity for spherical top- 
hat collapse at redshift z, extrapolated to today, Dl{z) is the linear growth factor, and cr^(Mmin) 
is the variance of the linear density field extrapolated to today, and smoothed on a mass scale 
-A^min- We treat C, and Mmin as free parameters. The parameter C, is an efficiency parameter that 
depends physically on the number of ionizing photons produced per stellar baryon, the fraction of 
gas converted into stars, the fraction of ionizing photons that escape from host halos and make it 
into the IGM, as well as the average number of recombinations during reionization. This expression 
is similar to the form Xe{z) = C/coiK-z) (with fcow denoting the collapse fraction in halos above Mmin) 
that is sometimes used, but our form goes smoothly to Xe = 1 as C/coU becomes large. The parameter 
Mmin is the minimum mass for a halo to host an ionizing source. In reality C and Mmin likely vary 
with redshift and environment, but we neglect this here in order to adopt a simple two parameter 
model. 
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Figure 13: Same as Fig. [T2l but for physical model with the minimum mass of Mmin = 2 x 10^ Mq 
and efficiency of C = 50. While CMBPol cannot distinguish between different physical 
models along the direction of having the same r value, it can distinguish between the 
physical model and an instantaneous reionization model with the same r at ~ 95% CL. 



We computed a physical model which is tuned to have the same optical depth r = 0.12 as 
the high r model in the previous section; it turns out to have a roughly comparable spread in 
the reionization epoch as well. Fig. [13] shows that CMBPol can distinguish between the physical 
model and an instantaneous reionization model with the same r at ~ 95% CL. As in the Zr — A 
parametrization, it is hard to distinguish models that have the same r, yet different durations: only 
the instantaneous model can be rejected. This is seen in the contour plot as a degeneracy between 
very efficient yet rare source models (high high Mmin) where reionization occurs rapidly, and low 
efficiency abundant source models (low ^, low Mmin) in which reionization is gradual. This supports 
the basic conclusion of the previous section: CMBPol can reject very rapid reionization models, but 
detailed constraints on the duration of reionization are impossible for smoothly varying reionization 
models. 

3.1.5 Constraints on Principal components 

The true reionization history may be more complex than the parametrized models of the previous 
two sub-sections. For example, feedback effects from the first generation of sources will likely impact 
the formation of subsequent galaxies and/or quasars, forcing the efficiency parameter and minimum 
host mass to vary with redshift and environment. If this is the case, one can hope to learn more 
about the reionization history than implied by the previous two sub-sections. In addition to the 
timing of reionization, and constraints on the high redshift contribution to the opacity, one may 
hope to learn further information regarding the duration of the reionization process. For instance, 
is the reionization process even monotonic? 
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Figure 14: Left: Ionization histories for the fiducial model used for MCMC with principal compo- 
nents (solid) and for an instantaneous reionization model with the same total optical 
depth, T = 0.077 (dashed). Right: E'-mode polarization power spectra for the ionization 
histories in the left panel. The 68% CL errors for a cosmic variance limited CMBPol 
experiment are plotted for the fiducial model (binned in i). 



As an example of a more general reionization history, we consider a double-peaked reionization 
history from Furlanetto & Loeb [81[. These authors investigate the conditions under which double- 
peaked reionization histories are physically plausible, and find that such histories require fine-tuning 
and are somewhat difficult to arrange, but not impossible. We hence adopt their model with a 
clumping factor of C = 3 and a high virial temperature due to photoionization heating of Th = 
2.5 X 10^ K. The ionization history in this model, its ii^-mode power spectrum with CMBPol error 
estimates, and the same for an instantaneous reionization model with the same r (r = 0.077) are 
shown in Figure [TH The instantaneous and double-peaked model can be distinguished at high 
significance (in excess of 3-0"). This illustrates that, if more freedom is allowed in the underlying 
reionization history than in the previous sub-sections, one can constrain more than just the total 
optical depth. 

This motivates a more conservative approach that complements the constraints from the previous 
sub-sections. Here we allow Xe(z) to be a free function of redshift and see what constraints the data 
place on the form of this function with minimal theoretical assumptions. One simple implementation 
of this approach is to parametrize the ionization history using the values of Xg in several wide redshift 
bins 76|, l82|. An alternative parametrization that we employ here is principal components (PCs) of 
the ionization history, an orthogonal set of basis functions for Xe(z) ranked in order of how well they 
can be measured with large scale polarization data [83]. The interpretation of constraints on PCs 
is less intuitive than for the binned ionization fraction, but the advantages of PCs are that they are 
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largely decorrelated from one another, there is less need for arbitrary redshift bin choices, and the 
impact of reionization on the CMB polarization is concentrated within a small number of the best 
measured components. 

The principal components S^{z) are defined as the eigenfunctions of the Fisher matrix that 
describes the information about reionization contained in the CMB polarization power spectrum on 



large scales [83|, 



^ V 2y dXeiZi) dXe{Zj) 

where ^max ~ 50 — 100 is large enough to include all effects of the mean ionization fraction, and the 
factor [Nz + 1)"'^ is included so that the eigenfunction shapes are independent of the initial redshift 
bin width Az = Zj+i — Zi as Az 0. The ordering of the PCs is based on their expected variances, 
given by the inverse eigenvalues of Fij, cr^, with the best constrained PCs having the lowest values 
of a^. For the set of PCs used here, we compute Fij at a fiducial model of constant Xg*^ = 0.15 
in redshift bins with spacing Az = 0.25 over a redshift range 6 < z < 30. The first five of these 
components are plotted on the left side of Fig. [T5j The = 6 is motivated by 

observations of quasar absorption spectra. For details on how results depend on the chosen fiducial 
model and maximum redshift, z^g_^, see Ref. [To']. The above equation assumes full sky coverage and 
sensitivity limited only by cosmic variance, but these choices have little effect on the PC shapes. 
Any function of redshift between Zjam and Zmax can be expressed as a linear combination of PCs, 

N 

Xe{z)=xf{z) + Y,m^^S^{z). (5) 
Only the 3 — 5 lowest variance PCs affect the polarization power spectrum to a significant degree 



(i.e. larger than or comparable to cosmic variance) [79|, |83|]. One can therefore truncate the sum 
in Eq. ([5]) when estimating parameters from CMB data with = 3 — 5 while still retaining the 
full information about reionization present in the large scale polarization. For the results presented 
here and in § 13.1.21 we always use the five components plotted in Fig. [15] as parameters for MCMC 
analysis. 

The right panel of Fig.[T5]shows MCMC constraints on five PCs and the total optical depth r from 
current and future CMB data. The largest set of contours are from five-year WMAP temperature and 
polarization power spectra, using MCMC estimates of the five reionization parameters in addition 



to the baryon density, CDM density, Hubble constant, and scalar amplitude and tilt 75|. Forecasted 



constraints on the five PCs from polarization at ^ < 50 are plotted for Planck, assuming a sensitivity 
of Wp^^'^ = 63.1 //K arcmin and beam size 0fwhm = 7 arcmin, and for a high resolution, cosmic 
variance limited experiment ("CMBPol"). The fraction of sky usable for cosmological constraints is 
assumed to be /sky = 0.8 for both future experiments. We include the temperature data for WMAP 
since it provides some additional constraints to the current polarization data, but for Planck and 
CMBPol essentially all of the information about the mean ionization history can be obtained from the 



large scale polarization [75l. l79l|. To find the projected Planck and CMBPol parameter uncertainties 
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Figure 15: Left: Five lowest variance principal components of the ionization history for large scale 
£^-mode polarization. Right: Constraints on principal components and total optical 
depth to reionization from 5-year WMAP data including the full temperature and 
polarization power spectra {blue), and forecasts for Planck {gold) and a cosmic variance 
limited CMBPol experiment {red) using polarization data at ^ < 50. Each set of 
contours shows 68% and 95% CL regions. 

we include polarization data up to ^ = 50 since treating the ionization history as a free function up to 
high redshift allows for significant variation in Cf^ at smaller scales than in typical parametrizations 
of reionization. For these forecasts, the amplitude of the temperature and polarization spectra at 
high i is held fixed to the best fit to current data by keeping the value of A^e"^"^ fixed for each 
MCMC sample. 

While the WMAP constraints in Fig. [15] come from real data and therefore reflect the true 
reionization history, we must choose the model for the simulated data used in the Planck and 
CMBPol forecasts. We take the fiducial ionization history to be the double-reionization model of 
Figure [T3] and we compute MCMC likelihoods relative to the E-mode polarization spectrum of this 
model. This fiducial model has a total optical depth of r = 0.077 and a period of reionization 
extending from z ~ 20 to z « 7. The estimates of r and the principal component amplitudes 
recovered from the MCMC analysis agree well with the true parameter values of this input model. 

The boundaries of plotted regions for the principal components on the right side of Fig. [15] are the 
edges of our priors on individual PCs. All potentially physical ionization histories (i.e. < Xg < 1) 
are included within these priors, including those that appear to violate the allowed range of Xg when 
the sum of PCs is truncated [tO^]. However, this does not mean that every model inside these bounds 
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is physically allowed or plausible. The PC priors are therefore in the spirit of the conservative 
approach to reionization constraints in which we allow for the largest conceivable variation in the 
ionization history. 

If we allow the priors on PC amplitudes to define the range over which measurements enable 
interesting comparisons between models, then we find that current data places tight constraints on 
the first principal component, roughly corresponding to r, and weakly limits the allowed range of 
the second component. Planck and CMBPol are expected to provide accurate estimates of several 
components: about 3 — 4 for Planck and 4—5 for cosmic variance limited CMBPol. Improvements 
in T and the two lowest variance components between Planck and CMBPol are minimal, with most 
of the impact of having truly cosmic variance limited data appearing in the higher PCs. 

It is important to note, however, that principal component parameters estimated at a significant 
level compared with the conservative < Xe < 1 priors do not necessarily lead to interesting 
constraints on physically motivated models. For many models of reionization that are currently 
thought to be realistic, the variation between models in higher variance PCs that undergo several 
oscillations in redshift is likely to be negligible. Therefore, the number of PCs that can be measured 
to an "interesting" level of precision depends on the range of possible variation in each of the 
components for the set of models being considered. Although a cosmic variance limited experiment 
is expected to measure 4 — 5 components at a significant level compared with the full possible range 
when the ionization fraction is allowed to be a free function of redshift, under more restricted classes 
of models like those in ^ 13.1.41 it may be difficult to measure even two parameters at high significance. 

Until observational limits on reionization and our theoretical understanding greatly improve our 
confidence in the specific behavior of the ionization fraction over time, it is worth keeping in mind 
that the true reionization scenario might be stranger than what we expect. Forecasts for models with 
a solid physical basis provide a realistic view of the advances we can expect from one experiment 
to the next, while the projected uncertainties on parameters like principal components that form a 
more model-independent description of the ionization history place an upper limit on the information 
about reionization that can be extracted from the CMB power spectra on large scales. 



4 Impact of improved knowledge of reionization on other 
parameters 

The T determination for Planck will not be significantly degenerate with any other parameter apart 
from the amplitude of scalar perturbations and in some instances with the scalar-to-tensor ratio r 



7a]. This is also true when a reconstruction of the reionization history from the data is attempted 



however discrimination between specific reionization scenarios may be a challenge for Planck |84i |. 

As for degeneracies, in such experiments r would only be significantly degenerate with the 
amplitude of scalar perturbation Ag and as a consequence with erg (see FigllGp. Even in non-minimal 
scenarios where an arbitrary scalar running spectral index and an free tensor spectral index are 
considered, r does not show significant degeneracy with these parameters (Fig llTp . 

If CMBPol is able to measure the i?-mode polarization signature of gravitational waves from 
inflation, accurate estimates of reionization parameters will be necessary to avoid biasing constraints 
on inflationary parameters. Unless the amplitude of tensor perturbations is close to the current upper 
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Figure 16: Joint confidence regions, at 68% and 95% CL, in the t-Aj and r-crg, assuming an 
instantaneous reionization. Curves are plotted for Planck (green lines) and CMBPol 
(red). There are no significant degeneracies between r and the remaining parameters. 




Figure 17: Same as Figure [T6] but for a fiducial model with a non-zero running of the scalar spectral 
index, rirun and a non-zero tensor spectral index, nt- 



limit set by the CMB temperature data, the ability to probe inflation is likely to rely heavily on 
the reionization peak in the large scale i?-mode tensor spectrum since instrumental noise, angular 
resolution, and the conversion of E'-modes into i?-modes via lensing all present significant challenges 
for measuring Cf^ on smaller scales. Without sufficient information about the reionization history 
from large scale E-mode polarization, estimates of infiationary parameters such as the tensor-to- 
scalar ratio r and tensor tilt nt that make use of the reionization peak in Cf^ are strongly correlated 
with the total optical depth to reionization. Accurate constraints on r and on the shape of the 
reionization history from cosmic variance limited measurements of Cf^ would break this degeneracy, 
enabling the use of the large scale tensor contribution to Cf^ for testing models of infiation 



5 Probing reionization with secondary anisotropies 
5.1 Small angular scale measurements 

On scales comparable to a tenth of a degree or smaller (i > 2000), a second reionization observable 
comes into play, anisotropy created by the inhomogeneous distribution of ionized hydrogen (HII) 
surrounding the first sources. In this section we will discuss this probe and its potential to constrain 
not only the mean reionization redshift, but also its duration. Whether CMBPol or future ground 
based experiments will be the best platforms to do this measurement is currently unclear. A po- 
larization mission which attempts to de-lens the polarization signal to reduce the lensing B-mode 
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contamination requires high angular resolution and high sensitivity. Such a mission, if also sensitive 
to temperature anisotropics, would be ideal to study the secondary anisotropics generated during 
reionization. 

There are two observables sourced by the inhomogeneous distribution of ionized hydrogen during 
the EoR. First note that anisotropics due to varying integrated optical depth T(n) lead to anisotropics 
in the observed power in different directions on the sky. Because in typical reionization models the 
sightline fluctuations are only of order 1%, the contribution to the power spectrum is of order 10~^, 
which is below cosmic variance. 

The dominant small scale CMB reionization observable comes instead from the fact that the 
extended ionized bubbles (which are about 30Mpc/h across at the intermediate EoR stages, lead- 
ing to typical signatures of size 10' on the sky) are embedded in large scale line-o f-sig ht velocity 
flows. While the standard Ostriker-Vishniac/kinetic Sunyaev-Zel'dovich (SZ) effect |86l. l87l|. owing 
to homogeneously ionized gas, is suppressed at these high redshifts, due to relatively low gas temper- 
atures and small density contrasts, the inhomogeneity of ionized hydrogen adds a large contrast to 
the signal. The effect has been modeled in various ways [88|, |89|, |90|, l9l|, |92] • In [9l|, |92] the signal was 
found to peak at i = 2000, where it is still almost two orders of magnitude below the primary CMB 
signal. It is approximately one order of magnitude below the primary fluctuations at £ = 3000, and 
becomes the dominant contribution over the lensed primary CMB at i = 3500 — 4000, depending 
on the details of the reionization model and on how well the thermal Sunyaev-Zel'dovich effect [93 1 
can be removed from the data (see Figure [TH]) . The height of the "patchy reionization bump" (to 
be contrasted with the large scale bump in polarization discussed in the last section) depends on 
the total extent of the epoch, and its width depends on how much time the ionized fraction spends 
in various phases. If, for example, reionization stalls in its late stages, there will be relatively more 
large bubbles. 



In 



911] it was shown that already at Planck sensitivity patchy reionization can bias inferred 



cosmological parameter values if it is ignored, in particular if only temperature data are used. In 
polarization, secondary anisotropics are expected to be suppressed (see e.g. [89|, |9j]) by 4 orders 
of magnitude compared to the primary component. We make the approximation that the patchy 
reionization signal is unpolarized. The reason temperature data can be biased is that the broad 
bump seen in figure [18] enhances the power on small angular scales and, if not taken into account 
in the modeling, leads to (for example) an inferred primordial fluctuation spectrum that is bluer 
than it really is. The bias can be avoided either by only measuring temperature fluctuations out 
to ^ = 1500, or by including an extra parameter for reionization to marginalize over, hoping that 
this parameter will not degrade the constraints on other parameters much. We go the second route 
here, and investigate parameter constraints possible with full CMBPol temperature and polarization 
information. 

In this calculation we assume a raw detector noise level of Wrp = l/;xK-arcmin and ^fwhm = 3'. 
We model reionization following the fast semi-analytic scheme of 66|, [9l|. Here we applied this 
model in a 1 comoving Gpc/h simulation volume with 1,200^ volume elements. We also go beyond 
previous work in that we vary the threshold for overdensities to be ionized throughout the box while 
translating between comoving position and redshift, thereby creating a model of "reionization on 
the light cone". The redshift extent of the simulation box covers the entire reionization process. 
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Even if one attempted to measure this contribution in one broad band power. 
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1000 10000 

1 

Figure 18: Impact of patchy reionization on the small scale CMB anisotropics. The thick lines 
show contributions to the small scale CMB power spectrum assuming the spectrally 
distinguishable signal of thermal SZ and point sources can be "cleaned" . The patchy 
reionization excess over regular kSZ peaks at scales I ~ 2000. On scales smaller than I ~ 
4000, the primary CMB vanishes; however, cluster gas physics will have to be reliably 
modeled in order to discern the patchy reionization contribution from low redshift 
sources. 
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Modeling such a large volume in a full radiative transfer simulation would be prohibitive because of 
the requirement of a large dynamical range to model both the sources and large scales of reionization. 
However the size is particularly important for modeling the kinetic SZ effect because the velocity 
power spectrum peaks on scales of 250 Mpc/h. Our fiducial model, which will be used as template 
in the Fisher matrix analysis below, has a value for the ionization efficiency of 15, which leads in 
our concordance cosmology to an ionization fraction of 0.5 at z=ll. To obtain a template for the 
patchy kSZ signal, we integrate the kinetic SZ signal 



TCMB 

with 



;n) = „) y d>,|o-2e-'"<'''ic.('))lfl ■ 9 (6) 

where the Thomson scattering optical depth is ax, T^i is the optical depth from the observer to 
conformal time ry, and n is the line of sight unit vector. 

Assuming that the total duration of the EoR is long compared to the scale corresponding to 
the peak of the velocity power spectrum, we can make the approximation that if reionization takes 
e.g. twice as long to proceed, the r.m.s. of the induced anisotropy will be \/2 larger and the power 
spectrum amplitude doubled. Therefore to ffi'st order a 10% constraint on the amplitude of the 
patchy reionization admixture to the total CMB power spectrum can be translated into a 10% 
constraint on the duration of the process. 

We assume here that the thermal SZ component, which is roughly comparable to the kinetic 
component on large scales (see Fig. [TSl) . can be cleaned sufficiently. This can be justified by the fact 
that at ^ ~ 2000, which is where the patchy reionization contribution peaks, most of the thermal SZ 
signal stems from massive and hot clusters of mass M ~ 5 x lO^^Mpc/h and larger at low redshift 



(see e.g. Figures 5 and 6 of [95|). These large bright objects can hopefully be excised using the 
multi-frequency information of the experiment as well as X-ray data from ROSAT or planned 
missions like the Constellation-X or X-ray Evolving Universe Spectroscopy (XEUS) missions [^ . 
Furthermore, the thermal SZ vanishes at 220 GHz, leaving an unobstructed view on the patchy 
reionization contribution to an experiment observing in this frequency band. 

Table [2] shows constraints that can be obtained on 8 cosmological parameters including a running 
of the spectral index as and the patchy reionization amplitude. One can see that Planck will not 



be able to place interesting constraints on the reionization model. In [9l|] an extended reionization 
model (their model C) was investigated which was able to bias Planck constraints more strongly and 
might be borderline detectable. 

In the case of CMBPol we find most parameter biases to be large (first column), since they are all 
much more tightly constrained. When the patchy reionization "nuisance parameter" is included in 
the analysis, we see that CMBPol should be able to place sub-percent constraints on this amplitude 
parameter, see table [2j 



^^The brightness of polarized foregrounds will make observations in multiple frequency bands, roughly 
covering the range most interesting for the SZ, a necessity for a future CMB mission aiming to study the 
primordial E and B mode signals. 

-"^•^For details on planned X-ray observato- ries see http://constellation.gsfc.nasa.gov/ and 
http://www.rssd. esa.int/index.php?project=XEUS. 
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P [nd. valuej 


^(/o), ReiExp (Planck) 


RerExp, /max = 3000 


ReiExp, /max = 2000 


V r\ r\ r"! a'\ 

T-^i [0.084] 


2.70 (6.97) 


2.73 


2.78 


[0.742] 


0.30 (0.89) 


0.31 


0.41 


[0.1 093] 


0.41 (1.15) 


0.44 


0. 56 


a;;, [0.02255] 


r\ -1-1 //^ o \ 

0.11 (0.68) 


0.12 


0.23 


Us [0.963] 


0.17 (0.47) 


0.18 


0.21 


As [2.410-9] 


0.38 (1.22) 


0.48 


0.52 


«s[0] 


0.20 (0.60) 


0.23 


0.33 


Apr[1] 


0.42 (119.0) 


3.22 


33.5 



Table 2: Forecasted constraints, for a full sky Reference Experiment with noise level IfiK — arcmin 
(\/2//K-arcmin in polarization) and a 1' FWHM beam and Planck from combined tem- 
perature and polarization information on 8 cosmological parameters including the patchy 
reionization amplitude parameter Apr, which crudely parametrizes the duration of the 
EoR. Also included are the results for the reference experiment for small scale cutoffs, 
which could be due to limited telescope/beam size or point source confusion. An /max of 
3000(2000) would correspond to a 1.5(l)m telescope primary at 220 GHz. 



Here we have not discussed how well CMBPol might be able to discern reionization models with 



different source properties. For example feedback processes [9a] or recombinations [97] could lead 
to relatively larger /smaller typical bubble sizes, imprinting the characteristic bump in Figure [18] 
on larger /smaller scales (and changing the slope at high i, where the reionization signal comes to 
dominate over the primordial CMB). Furthermore, the width of the patchy reionization "bump" will 
depend on how much time is spent in each partial ionization stage. 

A number of ground based observatories with high angular resolution, in particular the South 
Pole Telescope (SPT0 and the Atacama Cosmology Telescope (ACT|^ are presently beginning 
observations. These experiments will primarily target the CMB anisotropics on the scale of a few 
arcminutes, to constrain the thermal SZ and the kinetic SZ effects from clusters and large scale 
structure at lower redshift as well as gravitational lensing. They also may be able to detect the 
patchy reionization admixture on these scales. 

5.2 The effect of metals 

As the first sources turn on an ionize the universe metals are also produced for the first time. These 
metals eventually enrich the inter-galactic medium (e.g. [s^). Here we concentrate on the resonant 
scattering signature of metal enrichment of the inter- galactic medium by the first stars on the CMB. 
CMB photons interacting with a species X at redshift resonant transition of resonant 



frequency Vrs, show the effect of the interaction if z^s = i^rs/i^abs^ 1 S2|- The resonant scattering 



of this species generates an optical depth to CMB photons [lOOl ]: Txiz) oc frsnx{z), where nx 



is the number density of the X species and frs is the absorption oscillator strength, tx modifies 



^''The website for SPT is http://spt.uchicago.edu/ 
^^http: / /physics. princeton.edu/act 
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temperature anisotropies: Ay = e '^^ Aj;^,,,.^ + Aj-^^^ . In the limit of low metal abundance rx ^ 1, 



we can retain only linear terms in tx- Thus as described in [99|], the change in the anisotropies 
(At — ^Torig) reduces to a blurring factor {—Tx^Torig) P^^^ ^ "^^^ term which is important only in 
the large scales. The large scale signal is less clean to interpret than the high i signal because at 
low multipoles generation and suppression terms tend to cancel out. The high d. signal in multipole 
space is 99(]: Aa^^ = —Tx^fm^^ ^-^^ proportional to the primary anisotropies, with a well 

defined frequency dependence: in the presence of only one resonant species at Zrs the effect will be 
non-zero only at frequency Uo = i^rs/(l + -^rs)- Thus the metals contribution to the power spectrum 
is proportional to the primordial CMB power spectrum with constant of proportionality given by 

In a given experiment, the lowest frequency channel corresponds to the highest redshift where 
the metal abundance is the lowest. Increments of metal abundances can be measured by studying 
the d ifference map and the difference power spectra between channels. Signal-to-noise considerations 



10ll | show that the main limiting factor is the accuracy of the cross-channel calibration. 

A purpose-built experiment with the properties of the EPIC 2m set up, with accurate control 
of systematics and cross-channel calibration might be sensitive to changes in 3-12% (2-0") of the 
Solar fraction of Oxygen abundance at 12 < z < 22 for reionizaton redshift Zre > 10 (OIII 84 /i m 
transition) and to changes in N abundance of 60% its solar value for 5.5 < z < 9 (2-o")(NII 205 iJ,m 
transition). 

Appendix: Likelihood Function 

For computing the constraints from the WMAP 5-year data, we shall evaluate the likelihood function 
directly in pixel space, following Appendix D of Page et al. (2007). 

For forecasting parameter constraints from Planck and CMBPol, we shall evaluate the likelihood 
functions in harmonic space. The form of the likelihood function for a theoretical model, Ci, given 
the data, Cf"^"-, is 

Imax 

-21nL(Q|Cf'^*«) = /..,5](2/ + l){ln[CrC7r-(Cm 

1=2 

fjEE lQTT,data . ]\[TT\ , fjTT f (jE E ,data j\^EE\ _ 2(jT E E ,data 

, I \ I I / / \ I I I I I / Q\ 



where Ni are the noisebias (which is assumed to be constant over multipoles) given by 



= I /^i^^str, (9) 




^F = It7^T:7^^^=1 ^K-^-ir. (10) 

— 1/2 —1/2 

Here, and Wp are the temperature and polarization sensitivities, respectively, which are 
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related to "AT/T per pixel" in Table 1.1 on pp.4 of the Planck Blue Book via 

^ 2.725 K X (AT/T per pixel for Stokes / in units of nK/K) 

X (Angular Resolution in arcmin), (11) 
u;-V2 ^ 2.725 K X (AT/T per pixel for Stokes Q & J7 in units of jiK/K) 

X (Angular Resolution in arcmin). (12) 

Note that we have ignored the effect of beam smearing, as we focus only on low multipoles where 
the beam effect may be ignored safely. 

In our calculation we integrate the likelihood function up to Imax = 30, beyond which little 
information is available for reionization. When we compute the likelihood function for CMBPol, we 
shall assume that CMBPol is limited only by the cosmic variance up to Imax-, i-e., 
up to Imax = 30. On the other hand when we compute the likelihood function for Planck, we shall 
combine the sensitivities in 70 GHz (LFI), 100 GHz (HFI), and 143 GHz (HFI), such that 

= _^ I _^ ^ =, (13) 

w-^l-" = (14) 

Prom the Planck Blue Book we find (if^TO^Gifz' '^7ioo GHz-> '^■tiis GHz) ~ (179.3, 68.1, 42.6) iJ,K arcmin, 
and (^p 7o^Gif2) ■w^p 100 GHz^'^piM GHz) ~ (255.6,109.0,81.3) jj^K arcmin. Therefore, the combined 
sensitivities for Planck are given by 

w^^^^ = 35.4 fj,K arcmin, (15) 

Wp^^"^ = 63.1 ijlK arcmin, (16) 

and thus we find 

Nf^ = 1.06 X 10-^ hK"^ str, (17) 

Ni^E = 3.37 X 10-^ hK"^ str. (18) 
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